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Computational complexity theory contains a corpus of theorems and conjectures regarding the 
time a Turing machine will need to solve certain types of problems as a function of the input size. 
Nature need not be a Turing machine and, thus, these theorems do not apply directly to it. But 
classical simulations of physical processes are programs running on Turing machines and, as such, are 
subject to them. In this work, computational complexity theory is applied to classical simulations 
of systems performing an adiabatic quantum computation (AQC), based on an annealed extension 
of the density matrix renormalization group (DMRG). We conjecture that the computational time 
required for those classical simulations is controlled solely by the maximal entanglement found 
during the process. Thus, lower bounds on the growth of entanglement with the system size can 
be provided. In some cases, quantum phase transitions can be predicted to take place in certain 
inhomogeneous systems. Concretely, physical conclusions are drawn from the assumption that the 
complexity classes P and NP differ. As a by-product, an alternative measure of entanglement 
is proposed which, via Chebyshev's inequality, allows to establish strict bounds on the required 
computational time. 
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I. INTRODUCTION 

Global optimization is one of the most important com- 
putational problems in science and technology. But be- 
yond its practical relevance, it is also of deep theoret- 
ical interest when viewed from the broader perspective 
of computational complexity theory [l], Problems 
are ranged into an intrincatc classification by theoretical 
computer scientists, and an impressive corpus of theo- 
rems and conjectures has been built to relate them, such 
as the Cook-Levin theorem 0, 0] which proves the exis- 
tence of NP-complete problems, or the conjecture that 
P^NP. 

All those complexity classes are denned with respect 
to an abstract computer, the Turing machine. Physical 
devices designed to solve a particular problem need not 
be subject to that restriction, i.e.: a NP-complete prob- 
lem might be solved in polynomial time by a physical 
device even if P^NP. The reason is that Nature need 
not be a Turing machine. Notwithstanding, simulations 
of physical processes on classical computers are bound 
by the previous hierarchy of classes, since they are (ap- 
proximately) Turing machines. If P^NP, any attempt 
to solve a NP-complete problem in polynomial time with 
a simulation of a physical process on a classical computer 
must fail. The reasons for the failure must be deducible 
from the simulation details, and insight about the under- 
lying physical process might be obtained. 

Quantum mechanics provides the most promising 
physical attempt to outperform classical computation, 
and among all the quantum computational techniques, 
we will focus on adiabatic quantum computation (AQC), 
also known as quantum annealing The possibility of 

using AQC to solve NP-complete problems in polynomial 
time is one of the most exciting problems in quantum 



computation, but it is not addressed in this work. We will 
focus on classical simulations of AQC built upon matrix 
product states (MPS) Q. Concretely, we will analyse a 
technique based on an adiabatic extension of the density 
matrix renormalization group (DMRG) , published 
as quantum wavef unction- annealing (QWA) [lfj. Along 
with it, we will put forward and discuss the DMRG- 
annealing conjecture, which states that the efficiency of 
the QWA simulations of AQC is controlled uniquely by 
the maximal entanglement found during the process. 

Being a classical computational technique, DMRG- 
based simulations of AQC can never solve NP-complete 
problems in polynomial time, unless P=NP. Accepting 
Py^NP and the adiabatic DMRG-conjccture to hold we 
can put lower bounds on the behaviour of entanglement 
during AQC processes. In some cases, it allows us to 
predict the existence of quantum phase transitions. 

Our work, thus, puts under a different light ideas that 
are known in the field. In 2003 Vidal showed that a quan- 
tum computation involving a fixed amount of entangle- 
ment could be efficiently simulated using a matrix prod- 
uct representation 0], thus showing that an exponential 
speed-up was only possible if the MPS bond dimension 
grows fast with the input size. But, assuming that the 
ground state (GS) of a certain hamiltonian can be de- 
scribed as a MPS of fixed dimension, how hard can it be 
to find it? In 2006, Eisert showed that this problem can 
be NP-complete [ll| . Indeed both results are not hard to 
reconcile within the DMRG-annealing conjecture frame- 
work, as we will show. In 2009 Hastings proved that AQC 
with fixed gap in ID would never achieve an exponential 
speed-up II 21. based on his first rigorous proof of an area 
law in ID [13|. In recent years a new field is being built, 
known as hamiltonian complexity, which considers the 
computational complexity of performing classical simu- 
lations of quantum systems (see [l4j for a recent review). 
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This paper is structured as follows. In section |TT] 
we review the basics of adiabatic quantum computation 
(AQC) (or quantum annealing) [fall- Section |TTT] de- 
tails the quantum wavefunction annealing technique, an 
adiabatic extension of DMRG, and the DMRG- annealing 
conjecture is formulated and discussed in section HVl As- 
suming this conjecture to hold, our main result, which 
is the physical implications of complexity theory, are ex- 
posed in scction|Vj The paper closes with the conclusions 
and further work. 



II. ADIABATIC QUANTUM COMPUTATION 

Since the seminal article of R.P. Feynman in 1982 [Til ], 
physicists have had an increasing interest in the simula- 
bility of quantum mechanics, which has grown into the 
field of hamiltonian complexity [l4l |. The difficulties re- 
side in the exponential growth of the dimension of the 
Hilbert space. Quantum computation was born with the 
idea of converting this handicap into an opportunity: per- 
haps clever exploitation of this exponential growth will 
allow us to achieve an exponential speed-up of classi- 
cal algorithms, maybe to solve NP-complete problems in 
polynomial time |lq |. This hope has not yet been either 
fulfilled or disproved, and we will not address it here. 

Among the quantum computational techniques pro- 
posed, we will focus on adiabatic quantum computation 
(AQC) [3, studied also under the name of quantum an- 
nealing [a, AQC was proved in 2004 to be universal 
in the sense that the results of any quantum computation 
can be simulated in polynomial time with an AQC fl8l ]. 

An AQC is implemented by a physical device which 
establishes an adiabatic route between two hamiltonians, 
Hq and Hi, such that the ground state (GS) of Hq is 
easy to obtain physically, and the GS of Hi provides the 
solution to some problem. The GS of a hamiltonian is 
difficult to achieve experimentally when the system is 
subject to ageing, i.e.: when the low energy spectrum 
is complex, as it happens for most disordered systems. 
The adiabatic theorem ensures that, if the process is slow 
enough and the gap never vanishes, the ground state of 
Hi will be obtained from that of Hq. 

As a relevant example throughout this work let us 
consider the (classical) spin-glass problem [l9j]- Given 
a graph Q of N spins and a set of arbitrary real coupling 
constants Jy- attached to each graph link, we define the 
(classical) spin-glass energy as 

E = ~ £ Jno-io-j (1) 

(id) 

Where the Oi are values in { — 1, +1} attached to each 
site. The (classical) spin-glass problem is to find the 
values for Ui which minimize the previous energy. 

If the graph is ID, the problem is trivially in P. If it 
is 2D, a non-trivial construction found by Barahona (20j 



also renders the problem polynomial. For higher dimen- 
sions, or for random graphs of fixed connectivity (>3), 
the problem is NP-complete [2l[ . Even a 3D graph com- 
posed of two flat layers is NP-complete [20| ■ 

The AQC strategy for the spin-glass problem sets the 
destination hamiltonian, Hi, to be a quantum counter- 
part of eq. (H|), promoting the {— 1, +1} values of <x; to 
spin-1/2 operators (22|: 

H = - 2J JijSfSj (2) 

(id) 

To obtain Hq, we add to Hi a source of quantum fluc- 
tuations, such as a uniform transverse magnetic field: 
—T'^2 i 3f. Thus, T is the the tunable parameter which 
connects both hamiltonians and, as a function of it, 
we obtain the random Ising model in a transverse field 
(RITF) : 

#( r ) = -£ JijSfSj-r^sf (3) 

(id) 

where we see that H(oo) = Hq and #(0) = Hi. 

Let |^(r)) denote the ground state of the previous 
system as a function of T, which is only degenerate for 
r = 0. If r — > oo, the ground state is found just by 
making all spins point in the A-direction: 

|*(oo)) = -L(|+) + |-)f" (4) 

In this state, all classical configurations take exactly 
the same probability, so we may say that it is abso- 
lutely disordered. For T —> + , on the other hand, the 
ground state provides the solution to the classical spin- 
glass problem. Thus, the AQC strategy is to take T ^ oo, 
decrease it until F = 0, and then read the solution. The 
adiabatic theorem can be applied if the process is slow 
enough, assuming that the gap never vanishes. 

The main difficulty during an AQC experiment is to en- 
sure adiabaticity. The probability of a jump to an excited 
state increases exponentially as the energy gap closes, as 
reflected by the Landau-Zener formula J22|. Thus, if the 
system undergoes a quantum phase transition and the 
energy gap vanishes, the velocity must be reduced in an 
appropriate way at that point, increasing the computa- 
tional time. 

It may be tempting to try to extract conjectures about 
the minimal gap along an AQC trajectory from the (clas- 
sical) complexity class of the problem at hand. But 
these inferences are not valid, since the precise nature 
of the relation between the quantum and the classical 
complexity classes is not straightforward. Recent results 
of Altschuler and coworkers [23[ cast doubts on the pos- 
sibility of solving NP-complete problems in polynomial 
time using quantum computation, due to the very nar- 
row gap distribution in disordered systems which can be 
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deduced by Anderson's theorem. Nonetheless, other au- 
thors are more optimistic, believing that a route which 
avoids exponentially small gaps is feasible [13]. Those 
problems, which are of uttermost importance for quan- 
tum computation, will not be studied in this work. 

A caveat is in order: AQC is not designed for cases 
in which the gap is exactly zero at some point. This 
imposes certain restrictions on Hq and H\. For exam- 
ple: an adiabatic calculation in which we evolve from the 
(classical) ferromagnetic Ising model in S z to a (classical) 
spin-glass, also in S z , is not possible. Both hamiltonians 
commute, the ground state becomes degenerate at several 
moments, and the adiabatic theorem does not apply. 

III. CLASSICAL SIMULATION OF AQC: 
QUANTUM WAVEFUNCTION ANNEALING 

In order to apply the results of computational complex- 
ity theory, we should analyse algorithms running on Tur- 
ing machines, not on arbitrary physical devices. There- 
fore, we will study simulations of adiabatic quantum com- 
putation running on a classical computer. 

A first simulation approach to AQC is the use of path 
integral Monte-Carlo methods (PIMC) [12]. This tech- 
nique does not suffer from Landau-Zener level crossings, 
and the closing of the energy gap doest not constitute a 
problem. Nonetheless, if an attempt is made to solve a 
NP-complete problem using it, it is always found that, 
at some moment, the system undergoes critical slowing 
down. This forces long relaxation times and reduces the 
efficiency of the procedure. The exact amount of this re- 
duction is not easy to quantify, due to the different com- 
plex ity classes associated with probabilistic computation 

A different simulation procedure, quantum wavefunc- 
tion annealing (QWA), is a fully deterministic classical 
algorithm and lends itself more easily to analysis [Ioj |. 
The key feature of QWA simulation is that it computes 
the full wave function of the involved ground states. Let 
H(s) = (l-s)H + sHx, with s e [0,1]. Then, the QWA 
procedure is: 

1. Let s = and find the GS of the initial hamiltonian 
H(s) = H , |tf(0)>. 

2. Increase the adiabatic parameter: s — > s + Ss. 

3. Find the GS of H(s + 5s), |*(s + 5s)}, using the 
previous ground state as a seed. 

4. If the overlap | (^f(s + <5s)|\P(s)) | is below a given 
threshold, halve 5s and return to 2. 

5. If s + 5s < 1, make s = s + Ss and return to 2. 

If this computation were done in a naive way, the num- 
ber of stored components would be 2 N , thus it is unfeasi- 
ble. Instead, the wavefunctions may be stored as matrix 
product states (MPS): 



1*)= £ Tr(A^---A s -)\ Sl ,--- ,s N ) (5) 

S1---SJV 

where the A Si are 2N matrices of size to x m, and to is 
called the bond dimension. The total number of compo- 
nents in a MPS is, therefore, 2Nm 2 . Of course, m must 
be chosen so that the ground state is always accurately 
represented. The required value of the bond dimension 
to will be, therefore, of uttermost importance in order to 
evaluate the efficiency of the procedure. 

Our technique of choice in order to determine the MPS 
representation of the ground state for each value of T is 
the density matrix renormalization group (DMRG) [H, |9| . 
The technique has several important features: 

• DMRG is variational within the MPS subspace of 
the Hilbert space. 

• DMRG allows for adaptable values of to. In QWA 
simulations, to must be an open parameter, which 
is chosen to be large enough to represent the state 
accurately at each simulation stage, to a given tol- 
erance. 

• DMRG benefits from the use of a seed state in order 
to accelerate convergence, via the so-called wave- 
function transformations [261 ] . 

Thus, the QWA algorithm under consideration is just 
an annealed extension of the DMRG. 



IV. THE DMRG- ANNEALING CONJECTURE 

We conjecture that QWA simulations will always ob- 
tain the true optimum of the problem in linear time in 
the size and polynomial time in the maximal bond di- 
mension. Of course, the maximal bond dimension will 
itself depend also on the system size. The statement 
about linear time in the system size assumes that the 
bond dimension is held constant. In this section we will 
discuss the arguments in favour of this conjecture, and 
its possible pitfalls. 

We should stress that the bond dimension to is an 
adaptable parameter of the QWA computation, which 
will always take the minimum value that makes the rep- 
resentation of the wavefunction accurate up to a certain 
tolerance. The maximal to is defined to be the highest 
value found during the full computation. This value will 
be, in general, highly dependent on the DMRG path, i.e.: 
the ordering of the spins required for the sweeping. Ob- 
viously, it is desirable to choose a path that minimizes 
this maximal to. This minimization problem can be also 
a hard problem, as we will discuss in section fVl Heuristic 
low-cost approaches have been attempted to find a good 
DMRG path, with reasonable success [27| ■ 

It is usually assumed in the DMRG literature that the 
bond dimension to scales as the exponential of the von 
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Neumann entropy of entanglement @, m w exp(S). It is 
straightforward to show that, given a MPS of dimension 
m, the maximal von Neumann entanglement that can be 
supported is indeed log(?n). But this scaling can only 
be proved rigorously for some mild distributions of the 
eigenvalues of the density matrices. 

The eigenvalues of a reduced density matrix p consti- 
tute a discrete probability distribution, {pi}^i , which we 
will assume to be in decreasing order. The approximation 
of assuming a certain bond dimension m is equivalent to 
retaining the first m eigenvalues pi and neglecting the 
rest. Given a certain tolerance e > 0, we would like to 
find m(e) to be the number of eigenvalues that must be 
retained so that the sum of the remaining ones is smaller 
than e: 

E Pi < 6 ( 6 ) 

i>m(e) 

The von Neumann entropy of entanglement is the 
Shannon entropy of the eigenvalues of p: S = (— log(pj)) . 
If the eigenvalues decay exponentially, it can be proved 
that for all e > 0, m(e) oc e~ 1 exp(S'). In case of slower 
decays, still a polynomial relation can be found between 
cxp(S') and m, but this is not true for generic distribu- 
tions. 

Instead of the von Neumann entropy, we can employ a 
different measure of entanglement, such as the variance 
of the eigenvalue index: 

n t / n t \ 2 

a* S £i». ft - 5>A ^ 2 )-« 2 (7) 

i=l \i=l / 

With this new measure of entanglement, Chebyshev's 
inequality directly provides the desired relation: 

m(e) = (i) + -^<r => V k<( (8) 

V 6 -IS 

v i>m(e) 

Therefore, whenever a is finite, we have a rigorous 
bound on the bond dimension as a function of the toler- 
ance. 

At each step, the QWA tries to transform the GS at a 
value of s into the GS at a slightly larger value, s + Ss. 
In order to perform this transformation, it will proceed 
by carrying out as many DMRG sweeps as necessary. 
The cost of each DMRG sweep scales linearly in N and 
polynomially in m. But there is still one caveat: the 
number of DMRG sweeps might scale too fast with the 
system size. The conjecture assumes that this will not 
be the case, since both states |\&(s)) and |\&(s + Ss)) will 
always be close enough. The fidelity | (^ r (s)|4 f (s + Ss)) \ 
measures the speed of change of the wavefunction [28| . 
Evidently, the algorithm must proceed more slowly near 
the quantum phase transition. 



Thus, the strongest claim within the DMRG- annealing 
conjecture is that the number of required sweeps will al- 
ways scale polynomially with the system size. Until a 
proof (or a counterexample) is provided, we can only re- 
fer to the community consensus that the DMRG is sur- 
prisingly robust in this regard p9| . 

To summarize: there are two steps in the QWA com- 
putation which depend on a certain tolerance threshold: 
the DMRG tolerance (on the sum of neglected probabili- 
ties) and the annealing tolerance (on the minimal allowed 
value for the fidelity). The DMRG tolerance is strictly 
controlled by entanglement, via the bond dimension. The 
effects of the annealing tolerance, on the other hand, have 
not been explored in depth, although no surprises are ex- 
pected. 

Some recent works have been devoted to study the ef- 
ficiency of the computation of MPS. Finding the ground 
state of a ID quantum hamiltonian can be even NP- 
complete [29[. If m is kept fixed, MPS are always nicely 
approximable [30| . This means that it is always possible 
to obtain an approximation, within a factor (1 + e) of the 
ground state energy, polynomially both in time and in 
e. The complexity class which conveys fact this is called 
full y polynomial-time approximation scheme (FPTAS) 

There are other algorithms to simulate a quantum 
computation based on the MPS representation, e.g. Q 
for standard quantum computation or (32j for AQC. The 
latter case uses real time simulation, which may lead to a 
new source of loss of adiabaticity. Nonetheless, they both 
point to similar relations between entanglement and the 
efficiency of the calculation. 

Since the efficiency of the procedure is controlled by 
entanglement, via the maximal-MPS dimension, let us re- 
view the known results regarding its behaviour. Gapped 
systems usually fulfill the area law [HI, [H| , which states 
that the von Neumann entanglement entropy between 
two parts of a system scales as the number of broken 
links among them. This result was proved rigorously in 
ID by Hastings [l3|. In some higher dimensional inho- 
mogeneous and/ or disordered systems it is known not 
to hold [H, [36(. Moving to critical points, von Neu- 
mann and Renyi entanglement entropies have proved to 
be very useful magnitudes in order to pinpoint quantum 
phase transitions (QPT) [2?]]. During any adiabatic pro- 
cess with finite N, QPT always presents a local maximum 
in S. It has also been conjectured that a non-analyticity 
of S may constitute a good indicator of a QPT [37[ . The 
ID case is again rather special. At criticality, S will scale 
with the system size. In some cases, the size-dependence 
of S can be found via conformal field theory (CFT) [38j . 
It has been shown that, for many critical ID problems, 
S(N) w a \og(N) + f3. Thus, assuming that the bond 
dimension m w cxp(S'), it will grow polynomially with 
the input size. 
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V. PHYSICAL IMPLICATIONS OF 
COMPLEXITY THEORY 

A. General principles 

Throughout this section, we will adopt the notation 
that a hamiltonian belongs to a complexity class if the 
problem of finding its ground state belongs to that com- 
plexity class. Of course, the DMRG- annealing conjecture 
will be assumed to hold. 

A complexity class which is simpler than P is LIN, 
i.e.: the class of problems that can be solved in linear 
time. Let us consider an AQC connecting two hamil- 
tonians Hq and Hi, such that Hq is in LIN but Hi is 
not. E.g., Hi may be in P but with computational time 
scaling faster than TV, perhaps only N\og(N). Then, 
the bond dimension must diverge at some point during 
the adiabatic route. Typically, this will imply that the 
von Neumann entropy will diverge also. This divergence 
is typically the hallmark of a quantum phase transition 
(QPT). 

The reason can be stated as follows. Let us assume 
that the entanglement entropy stays finite during the 
whole AQC procedure. Now, let us use a classical com- 
puter to run a QWA simulation of the AQC procedure, 
in time T(N). If the maximal bond dimension saturates 
with N, QWA results asymptotically in a linear algo- 
rithm to obtain the ground state of Hi, against the as- 
sumption. Therefore, entanglement must grow without 
bound with the system size at some moment during the 
AQC procedure, pointing to a QPT. 

In the same line, if Hq is P and Hi is NP-complcte, 
and P^NP, then any AQC connecting the two hamilto- 
nians will find, at some moment, a state with maximal 
bond dimension growing faster than polynomially in N. 
Typically this implies that the von Neumann entropy will 
grow faster than logarithmically with N. This state may 
correspond to a QPT. Again the reason is easy to state: 
otherwise, the classical simulation will provide a polyno- 
mial algorithm to solve a NP-complcte problem. 

In general terms, we may say that the adiabatic con- 
nection of two hamiltonians with different complexity 
classes puts restrictions on the physics along the path. 
In order to avoid violations of the results from complex- 
ity theory, entanglement must diverge at some moment 
during an AQC procedure. This divergence may take 
place as a quantum phase transition of a certain kind. 
It can be regarded as a kind of quantum censorship to 
prevent hard problems from being solved easily. 

This analysis is independent of whether we focus on 
average or worst-case complexity. Once the set of prob- 
lems is characterized, and a bound on the computational 
time is established, it can be immediately converted into 
a bound for the entanglement entropy for an AQC. 

A relevant point to be made is how to know whether 
the divergence of the maximal entanglement points to a 
quantum phase transition or not. In ID it is known that 
at a QPT both the entanglement entropy and the bond 



dimension diverge |13| . In higher dimensional studies, a 
more careful analysis is required, since the maximal bond 
dimension will likely diverge along the whole AQC pro- 
cedure. In any case, for any finite instance of a problem, 
the maximum attained by the entanglement determines 
the efficiency of the QWA. Therefore, if a QPT is present, 
even if bond dimension diverges in a full region, it is the 
scaling of entanglement at the QPT which establishes the 
running time for the simulation. 



B. Concrete examples 

Let us return to the spin-glass hamiltonian (J3J) - When 
r — > oo, the obtention of the ground state is a trivial 
problem, taking time 0(1). In ID, the obtention of the 
classical spin-glass minimum energy state is obviously in 
LIN. Therefore, our results do not apply in this case, 
since QWA takes always time > N . 

In 2D, on the other hand, a prediction can be done. 
Solving the 2D classical spin-glass problem is known 
to be in P [2(|. Therefore, the entanglement entropy 
must diverge for some value of T. We can only state 
that this maximal entropy will gro w, at least, logarith- 
mically. In fact, recent results [36( (cleverly exploiting 
the properties of the infinite randomness fixed point [39| , 
IRFP) show that it grows with a modified area law: for 
a block division cutting I links, the entropy scales as 
s(l) ~ Uog(log(Z)). Maximal entropy, as it is defined 
in this work, would be S(N) « N 1 / 2 log(log(/V)), thus 
rendering the time for the QWA simulation exponential. 
Our result is, therefore, valid but too weak. 

Nonetheless, the previous expression for the block en- 
tropy in a 2D quantum spin-glass is based on the average 
number of clusters cut by the block division. A well de- 
signed DMRG path might never cut more than one clus- 
ter at a time, just sweeping them one by one. In that case, 
the maximal entropy might grow much more slowly with 
the system size. But, in order to obtain such a path, one 
should first solve the classical problem. Therefore, again, 
our basic result is not violated. 

In 3D, or for random graphs of fixed connectivity, the 
NP-completeness of the problem forces the maximal en- 
tropy along the route to grow faster than log(TV). In this 
case, the result is not surprising. 

Other analysis [4(| have been carried out of the entan- 
glement entropy along typical standard quantum com- 
putations, and our general statements also hold. For 
example, the entropy along Grover's algorithm remains 
bounded, which is consistent with the fact that the prob- 
lem under consideration (unsortcd search) is in LIN. 
Also, an AQC designed to solve the exact cover prob- 
lem (which is NP-complete) found a QPT with S ps N. 
Shor's algorithm also shows a similar behaviour, although 
it is not clear which is the complexity class of the problem 
under study (i.e.: integer factorization). 

New predictions are easily made for AQC designed to 
solve problems which have never been studied. Thus, 
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an AQC designed to test planarity of a graph, or 2- 
colorability, need not find a quantum phase transition, 
since these problems belong to class LIN. But an AQC 
which sorts a set of numbers, or which performs the fast 
Fourier transform, will find a divergence in the bond di- 
mension and, very likely, in the von Neumann entropy, 
since their running time is larger than linear. The max- 
imal entropy in those cases might grow very slowly with 
size, since the (average) running time for the best al- 
gorithms are T « Nlog(N), so our only bound is that 
S should scale at least like log(log(A)). On the other 
hand, if P^NP, any AQC attempt to solve the traveling 
salesman problem, or 3-SAT, will always find maximal 
entropy growing faster than logarithmically. 

VI. CONCLUSIONS 

Theoretical physics has benefitted continuously from 
the incorporation of the results of branches of pure math- 
ematics which were born without any relation to it. Com- 
putational complexity theory is just another branch of 
mathematics, and this work just attempts to extract its 
most straightforward consequences for physics. As Niko- 
lai Lobachcvski put it: "There is no branch of mathemat- 
ics, however abstract, which may not someday be applied 
to the phenomena of the real world" [4l|. 

In this work we have put forward a strategy to derivate 
physical inferences from computational complexity the- 
ory. If a physical process is devised in order to solve 
some problem, simulating that process in a (classical) 
computer constitutes a (classical) algorithm to obtain 
the solution. The efficiency of this algorithm may be re- 
stricted by complexity theory, and this restriction must 
have some counterpart in the physical model which may 
apply to the real physical system. 

Concretely, adiabatic quantum computation (AQC) 
may be simulated in classical computers using quantum 
wavefunction annealing (QWA), which is a simulation 
strategy based on the density matrix rcnormalization 
group (DMRG). The efficiency of QWA is conjectured 
to be controlled by the maximal entanglement attained 
during the physical process. Different measures of en- 



tanglement arc discussed (bond dimension and von Neu- 
mann entropy) and a new one is introduced (variance of 
the eigenvalue index). Arguments are given in favour of 
this conjecture, along with an exposition of its possible 
pitfalls. 

If there is a bound on the scaling of the computational 
time to solve the problem on a classical computer, then 
this bound will transform itself into another bound for 
the maximal entanglement attained during the real phys- 
ical procedure. This type of no-go reasoning bears resem- 
blance to the second law of thcrmodyamics. In this way, 
the divergence of entanglement with the system size for 
some systems is proved. This divergence, in some cases, 
may be viewed as the apparition of a quantum phase 
transition, which can be regarded as a "quantum cen- 
sor" , preventing the solution of hard problems in an easy 
way (in a classical computer). 

The present derivation was performed using MPS and 
the DMRG, which are not specially well suited for mul- 
tidimensional systems, due to the need for a ID path 
to run through the system. Different generalizations of 
MPS exist, such as multiscale entanglement renormal- 
ization Ansatz [42[ (MERA) or projected entangled pair 
states [43[ (PEPS),which are altogether labeled as ten- 
sor product states [IJ] (TPS). New techniques have been 
developed for 2D optimization problems, making use of 
ideas related to dynamic programming [45L l46j . Also, 
other techniques have been proposed in order to simulate 
real time evolution in the Heisenberg picture We 
expect that application of this line of thought to these 
sophisticated tools will provide stronger predictions on 
the physics found during the performance of an adiabatic 
quantum computation. 
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